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Using the Gaussian model of the vortex tangle (VT) arising in the turbulent 
superfluid Hell, we calculate the energy spectrum E{k) of the 3D random 
velocity field induced by that VT. If the VT is assumed to be a purely frac- 
tal object with Haussdorf dimension Hp, the E(k) is a power-like function 
E{k) oc k~ 2+HD . A more realistic VT in Hell is a semi-fractal object, be- 
having as smooth line for small separations A£ <C R ( £ is the label coor- 
dinate, R is mean curvature )and having a random walk structure for large 
A£ with Hd = 2. For that case calculations give a spectrum E(k) that is 
k-independent for k smaller than 1/R (but larger than the inverse size of 
the system) and that scales as k^ 1 for larger k. The latter reflects the fact 
that for small scales a vortex filament behaves as a smooth line. Our results 
agree with recent numerical simulations. 
PACS numbers: 47.32.Cc, 47.37.+q, 67.40.Vs., 05.10.Gg. 



1. INTRODUCTION AND SCIENTIFIC BACKGROUND 

The idea that the properties of classical hydrodynamic turbulence can 
be described in terms of chaotically moving thin vortex tubes appeared some 
time ago (for a discussion and references see book by FrisclJil). Some addi- 
tional arguments in favor of that idea would be the recent direct numerical 
simulations of flowing liquids with the Reynolds number of order 100. It 
was shown that the vorticity field consists of a chaotic set of vortex tubes 
resembling the picture appearing in numerical simulations of the chaotic 
dynamics of quantized vortex filament in Hell!-! . Of course the visual 
similarity cannot be considered as a confirmation of the equivalence of both 
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pictures, therefore some quantitative analysis is required to check whether 
chaotically moving quantized vortex filaments do induce a 3D velocity field 
that possesses statistical properties similar to that of classical turbulence. 

In the paper we calculate the spectral density E{k) of the 3-D random 
velocity field induced by a chaotic set of quantized vortex loops. The spec- 
trum E(k) is one of the key characteristics of turbulent motion. To calculate 
it we use the Gaussian model of a VT elaborated by one of authors earlier.! 
We calculated spectral density E(k) for the cases of pure fractal lines with 
Haussdorf dimension Hp as well as for the more realistic case of a VT in 
superfluid turbulent Hell. 

2. GENERAL RELATIONS 



The average kinetic energy of the flow induced by a vortex loop can be 
evaluated as follows (seeO for details): 




Here s(£) describes the vortex line position parameterized by the arclength 
£, running from to the length of line L, s'(£) denotes the derivative with 
respect to arclength along the line (the tangent vector). Brackets () imply an 
averaging over all possible vortex loop configurations. Clearly, the integrand 
within the brackets () in (jlj) is just the distribution of energy dE/dk = 
Air^dE/d^k. in a spherically normalized k -space. 

Of course relation (||) is just a mathematical identity - the real physical 
meaning is hidden behind the averaging procedure. To fulfil an averaging 
procedure one needs knowledge of the statistics of a vortex loop. The latter 
should be taken from the appropriate stochastic theory of chaotic vortex 
filaments, which is unknown at present. Another way is to use numerical 
or /and experimental data. The latter gave rise to the Gaussian model of 
VT.U In that theory the statistics of a vortex filament is supposed to be a 
Gaussian. That supposition is quite natural if one realizes that a real vortex 
loop is formed as a result of numerous reconnections. Therefore the remote 
vortex line elements are weakly coupled. The only correlation between them 
appears due to propagating Kelvin waves. Then a step-by-step traversal 
along a line is a nearly Markovian process. Correspondingly, statistics of 
such filaments should be Gaussian with a sharply decreasing correlation 
function iV(£i — £2) = ( s '(£i) s '(£2)) between two remote tangent vectors. 

For Gaussian statistics the probability for the line to have some partic- 
ular configuration {s(£)} is expressed by the path integral containing the ex- 
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ponent of the bilinear combination of functions s(£). In practice, to calculate 
various averages it is convenient to work with the characteristic (generating) 
functional (CF), which is defined as a following average: 

W({P(0}) = /exp I • / P(£y(£) d£ I ) . (2! 



II 



Since the probability functional is a Gaussian one, calculation of the CF 
can be readily made by the full square procedure to give the result (here we 
consider only the isotropic case) 



L L 



W({P(0}) = exp - j J cWPte'MC " £")P(0 • (3) 







From the definition (§) of the CF it follows that the function N(£' - £") 
coincides with the correlation function (s^(£i)s^(£j)) between tangent vec- 
tors. Combining ([[J) and (||) onejjxmcludes that the spectral density E(k) is 
expressed via CF as follows (seeO for details): 



L L 2 

Here #(£) is a unit-step function. We will use relation (||) to calculate the 
spectral density of energy for various N(£' — £"). 

Before doing this let us discuss some fractal properties of vortex loops. 
Let us calculate the quantity D 2 = ((s a (£j) — s a (0)) 2 ), which is an averaged 
squared distance between initial (arbitrary) point of the curve s(0) and the 
points s(£j). Note that we deal with the real distance in three dimensional 
space (not along the curve!), therefore this consideration concerns the real 
size of the vortex loop embedded in 3D space. Using (||), the quantity 
((s a (£j) — s a (0)) 2 ) can be evaluated from relations: 

D 2 = dCi ^2<sL(£lK(6)> = / ?J f J ^1^2^(6-6)- (5) 

Jo Jo Jo Jo 

Let us assume that iV(£i — £2) is a power-law function 7V(£l — £2) = 
N * (£1 — £2)^- Here N is some constant with the proper dimension. Then 
the integral in the rhs of (||) is evaluated as oc £j +2 . That implies that the 
length L of the curve increases with its 3-D size D as L oc D 2 ^ x+2 \ The 
latter implies that the average loop is a fractal object having the Hausdorf 



(4) 
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dimension (HD) Hp = 2/(A + 2). For the example of a pure random walk 
chain with the "effective step" Iq, the function iV(£i — £2) = —£2), and 

the quantity A are equal (formally) to —1, Hd = 2, and the length grows as 
L oc D 2 . For a set of smooth lines A = 0, and Hd = 1. Vortex loops in the 
superfluid turbulent Hell are also described by Gaussian statistics. However 
the function iV(£' — £") is not a power-law function (seeH). The small parts 
of the filament behave as a smooth line whose (3-D) lengths are exactly equal 
to distance £' — £" along the curve. At the larger distance the filament has 
a random walk structure with an "effective step" equal to the mean radius 
of curvature R. A good approximation for the quantity N(£' — £") could be 
a function of type 1/(1 + (A£/R) 2 ). 



3. ENERGY SPECTRUM OF THE VELOCITY FIELD 

Let us calculate first the energy spectrum for the vortex filaments, whose 
statistics are Gaussian, with the CF ([|) and the function N(£' — £"), which 
is a pure power-law function JV(£' — £") = N * A£ A . Substituting this into 
d|]) one obtains an extremely cumbersome relation which in general should 
be evaluated numerically. However, since we intend to establish the connec- 
tion between the fractal properties of the vortex filament and the power-law 
spectra of 3D velocity field, it is possible to obtain some analytical results 
by focusing on counting powers. By use of the substitution z 2 = k 2 LA£ x+2 
relation (Q) reduces to an expression having the following structure. 

A+2 

kL 2 . 4+4A 2A \ 

E(k) = L J J Ak 2 (0 A+2 +B (j) X+2 \exp(-z 2 )z^k^dz (6) 

Here A and B are some numerical factors, of order of unity, appearing from 
calculating of integrals and E(k) is reduced by factor l/2irp s K 2 . It is seen 
that there are two different regions of the wave vector k, namely, large and 
small with respect to 1/L^ X+2 ^ 2 , which is nothing but the inverse 3D size 
of the loop. In the region k -C 1/ L^ +2 ^ 2 , the exponent is approximately 
unity, the term with B disappears due to the closeness of the line, and sim- 
ple counting of powers shows that E(k) oc k 2 . That distribution of the 
energy (valid far from the loop) describes an equipartition law and implies 
that dE/cfik is constant. In the opposite case of large k we can take the 
upper limit of the integral to be equal to infinity, after that integral is eval- 
uated as a combination of gamma functions (of constant arguments) mul- 
tiplied by /c - ( 2 + 2A )/( A + 2 ). That implies that the spectral density of energy 
E(k) oc k~ 2+HD . That result was obtained earlier from very qualitative 
considerations and is discussed in the book by Frisch.EJ 
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In parallel to a 'counting of powers' we performed a numerical calcula- 
tion of the spectra with the help of Mathematica 3.0. In Fig. la we depict (in 
logarithmic scale) three curves corresponding to the spectral density E(k) 
for a pure fractal vortex filament of the length L = 100 (all units are arbi- 
trary). The upper curve (see the region of small k) corresponds to the case 
A = 0, i.e. to smooth lines. It can be seen how dependence E(k) oc k 2 is 
changed to a dependence E(k) oc k~ l in the region of k ~ 100 -1 . That is 
a general result, since Hd = 1 and the 3D size of the curve is of order of 
100. The next (middle) curve corresponds to the case A = —4/5, i.e. to the 
so called self-avoid line. It is seen again that a £; 2 -dependence is changed 
to a dependence E(k) oc A; -1 / 3 in the region of k ~ 100~ 3/5 . Taking into 
account that Hd = 2/(A + 2) = 5/3, this behavior agrees with the general 
result. The lower curve corresponds to the pure random walk chain with 
j\T(£i - f 2 ) = W(6 - 6). Formally A = -1, and H D = 2. Therefore the 
& 2 -dependance is changed to /c-independent spectral density E(k) oc k in 
the region of k ~ lOO" 1 / 2 , as it should be. 

Let us now consider a VT in Hell. As has been discussed earlier, the 
VT in Hell is a semi- fractal object, behaving as a smooth line for small 
separation <C R, and having a random walk structure for large A£ with 
Hp = 2 and with the effective step equal to R. Therefore the speculations 
based on "counting of powers" are not applicable directly. However one can 
assert that the equipartition distribution of energy E(k) oc k 2 in the region 
of k « 1/v 7 LR should take place for a vortex loop in turbulent Hell. The 
region of large wave numbers is divided with some additional value k ~ 1/R. 
One can suppose that at that point behavior of the spectrum changes, ac- 
quiring the properties of pure fractal with the corresponding HD in the 
according regions. Namely, E(k) should scale as k in region 1/R << k, 
and be approximately independent l/V LR << k << 1/R (provided the 
latter is broad enough). Numerical calculations confirm this conjecture. In 
Fig. lb there is depicted (in logarithmic scale) the curves corresponding to 
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the spectral density E(k) for a vortex loop of length L = 100 with R = 0.01 
and R = 1. It is clearly seen that the upper curve (i? = 0.01 ) has three 
regions with bends at the points k\ = 1/V LR = 1, and k r = 1/R = 100. In 
three different parts we have E(k) oc k 2 , E(k) oc k° and E{k) oc k , corre- 
spondingly. That is in excellent agreement with qualitative considerations. 
When the mean radius of curvature is equal R = 10, the correspondent 
k\ = 1/VLR ~ 0.03, and k r = 1/R = 0.1. This behavior agrees with recent 
numerical simulationsQand our own. 

4. Conclusion 

Using the Gaussian approach for the chaotic vortex filament we calcu- 
late the energy spectrum E{k) of the 3D random velocity field, induced by 
that VT. If the VT is assumed to be a purely fractal object, the E(k) is a 
power-law function E(k) oc k~ 2+H °. The spectrum of a more realistic VT 
also agrees with that formula, if one thinks of VT as a semi-fractal object 
with a varying (depending on scale) HD. Let us note that the HD of a simple 
noncorrelated vortex loops satisfies the condition 1 < Hp < 3. That im- 
plies that Gaussian lines cannot create 3D flow, which has the Kolmogorov 
spectrum E(k) oc /c -5 / 3 . That implies that vortex loops should be strongly 
correlated (which, in turn, implies that the local approximation used in nu- 
merical calculations is invalid). Another variant is that the vortex filaments 
probably play no special role in turbulent flow. 
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